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1 Introduction 

Lattice computations have provided evidence for the existence of a phase transition 
in QCD [U [21 [3j H] notwithstanding the sign problem casts some doubt about [5]. 
From a theoretical standpoint, a missing low-energy limit obtained from QCD forced 
people to use some relevant models to cope with this question as the Nambu-Jona- 
Lasinio model or the sigma model. Starting from these, some authors proposed the 
introduction of an imaginary chemical potential [HI El IE]- Statistical confinement was 
proved in this framework and account was given for the phase diagram of QCD [9lfTU]. 
Recently, a different approach to identify a critical endpoint with a chiral chemical 
potential was also proposed [TTj . 

The aim of this contribution is to show how, provided a recent analytical approach 
to manage low-energy QCD, it is possible to give a proof of existence of a critical point 
for chiral symmetry in a limit at zero quark masses and chemical potential. The idea 
is to fix the form factor in a non-local Nambu-Jona-Lasinio model as discussed in [12] 
by obtaining it directly from QCD and deriving in this way the proper low-energy 
limit of the theory. Finally, we will be able to obtain a critical temperature in close 
agreement with lattice data. 

Our strategy will be to start from an analysis of classical field theory in the limit 
of a strong coupling and then, with the knowledge acquired in this way, we will be 
able to formulate a quantum field theory in this same limit. Extending the analysis 
to finite temperature will be accomplished with standard techniques. We will closely 
follow Ref.[13l. 



2 Classical field theory 

A classical field theory for a massless scalar field is given by 

d 2 <p + \<f) 3 = j. 
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The homogeneous equation can be solved exactly by 

<f> = fjL f^\* sa{p-x + e,i) (2) 

being sn an elliptic Jacobi function and \x and 9 two constant. This solution holds 
provided the following dispersion relation holds 

P 2 = ^\jl (3) 

so this solution represents a free massive solution notwithstanding we started from a 
massless theory. Mass arises from the nonlinearities when A is taken to be finite rather 
than going to zero as normally done in weak perturbation theory. In this latter case, 
it would be very difficult to recover the above solutions. When there is a current we 
ask for a solution in the limit A — > oo as our aim is to understand a strong coupling 
limit. So, we check a solution 

being 6<j) all higher order corrections. One can prove that this is indeed so provided 
50 = k 2 A J d A x'd A x"G(x - x')[G{x' - x")] 3 j{x') + 0{j{xf) (5) 

with the identification k — the same of the exact solution, and d 2 G(x — x') + 
X[G(x — x')] 3 = /i _1 5 4 (a; — x'). This implies that the corresponding quantum field 
theory, in a very strong coupling limit, takes a Gaussian form and is trivial (triviality 
of the scalar field theory in the infrared limit). All we need now is to find the exact 
form of the propagator G(x — x') and we have completely solved the classical theory 
for the scalar field in a strong coupling limit. In order to solve the equation 

d 2 G(x - x') + X[G(x - x')] 3 = ^-^(x - x') (6) 

we can start from d = 1 + case d 2 G (t-t') + \[G (t-t')} 3 = /j, 2 5(t-t'). The Fourier 
transformed solution is straightforwardly obtained as 

71-2 (_l\n e -(n+±)TT j 

GoM = 5> + l) K2[i) 1 + e _ (2w+1)ir U 2_ ml + i€ ( 7 ) 



n=0 



being m n = (2n + (f) 4 H and K(i) w 1.3111028777 an elliptic integral. We 

are able to recover the full covariant propagator by boosting from the rest reference 
frame obtaining finally 



oo 



7T 2 (-l)™ e -("+i^ 1 

2n+ lj- 

n=0 



G(p) = J2(2n + p2 _ < + te (8) 
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This shows that our solution given above indeed represents a strong coupling expan- 
sion being meaningful for A — > oo. 

Similarly, for a classical Yang-Mills field one has to solve the set of nonlinear 
equations 

d^Al - (l - ~) d^Aft + gf abc A b »{d,A c u - cU£) + gf abc d^(A b fM Al) (9) 

+g 2 f abc f cde A b ^Al = 

For the homogeneous equation, we want to study it in the formal limit g — > oo. We 
note that a class of exact solutions exists if we take the potential A® just depending 
on time, after a proper selection of the components (we dub this Smilga's choice [T¥j). 
These solutions are the same of the scalar field when spatial coordinates are set to 
zero (rest frame). Differently from the scalar field, we cannot just boost away these 
solutions to get a general solution to Yang-Mills equations due to gauge symmetry. 
Anyhow, one can prove that the mapping persists but is just approximate in the limit 
of a very large coupling. This mapping would imply that we will have at our disposal 
a starting solution to build a quantum field theory for a strongly coupled Yang-Mills 
field. This solution displays a mass gap already at a classical level. Exactly as in the 
case of the scalar field we assume the following solution to our field equations 

A a , = nj d"x'Df v {x - x')j hv {x') + 6A; (10) 

Also in this case, apart from a possible correction, this boils down to an expansion 
in powers of the currents as already guessed in the '80 [15J. This implies that the 
corresponding quantum theory, in a very strong coupling limit, takes a Gaussian form 
and is trivial. The crucial point, as already pointed out in the eighties |T5l [16], is the 
exact determination of the gluon propagator in the low-energy limit. This is possible 
using the following mapping theorem: 

MAPPING THEOREM: An extremum of the action 



S = Id x 



\m 2 - -/ 



is also an extremum of the SU(N) Yang-Mills Lagrangian when one properly chooses 
A ^ with some components being zero and all others being equal, and A = Ng 2 , being g 
the coupling constant of the Yang-Mills field, when only time dependence is retained. 
In the most general case the following mapping holds 



A;(x)=^(x) + 0(l/ViV 3 ) 1 (12) 
being rf some constants properly chosen, that becomes exact for the Lorenz gauge. 
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This theorem was proved in [T7] and, after considering a criticism by Terry Tao, 
in [18]. Tao agreed with this latest proof [T9] . 

The mapping theorem permits us to write down immediately the propagator for 
the Yang-Mills equations in the Landau gauge for SU(N): 



D ^ s - s^f^J?. +0 , , (13) 

being 

R = On 4-1] 



71-2 (_l)n+l e -(n+±)n 

B n =(2n + 1) T ^ \ (14) 



m„ = (2n + l) — — A. (15) 



and 

-m,„ = f9n -I- 1 1 

2K(i) 

The constant A must be the same constant that appears in the ultraviolet limit by 
dimensional transmutation, here arises as an integration constant [20] • This is the 
propagator of a massive field theory and the mass poles arise dynamically from the 
non-linearities in the equations of motion. But, we are working classically yet and we 
cannot claim anything about the quantum theory at this stage. 



3 Quantum field theory 



We can formulate a quantum field theory for the scalar field starting from the gener- 
ating functional 



z \j\ = / [#]exp 



(16) 



We can rescale the space-time variable as x — > \/\x and rewrite the functional as 

1 



exp 



A 



d*x 



(17) 



Then we can seek for a solution series as 



and rescale the current 



j — >■ j/X being this arbitrary. It is not difficult to see that the leading order correction 
can be computed solving the classical equation 



+ 



j 



that we already know how to manage. This is completely consistent with our preced- 
ing formulation [21] but now all is fully covariant. We are just using our ability to 
solve the classical theory. Using the approximation holding at strong coupling 



b = fi J d 4 xG(x — x')j{x') 



(19) 
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it is not difficult to write the generating functional at the leading order in a Gaussian 
form 



Z [j] = exp 



/ d 4 x'd 4 x"j{x')G{x' - x")j{x") 



(20) 



This conclusion is really important: It says that the scalar field theory in d=3+l is 
trivial in the infrared limit. This functional describes a set of free particles with a 
mass spectrum 



m n = (2n + 1) 



7T 



2K{i) 



AV 



2' " 



(21) 



that are the poles of the propagator, the one of the classical theory. We note that 
this propagator is describing free massive particles with a superimposed spectrum of 
a harmonic oscillator as they would have a structure. 

For Yang-Mills theory, the generating functional can be written down with the 
following terms in the action 



Sym = -\J d 4 xTrF 2 - -L J d 4 x(d ■ A) 2 , 



the ghost field 



(22) 



(23) 



S g = - J d 4 x(c a d fi d fl c a + g^f abc d^c c ) 
and the corresponding current terms 

S c = J d 4 xj^(x)A a »(x) + J d 4 x [c a {x)e a {x) + e a (x)c a (x)] . (24) 

We now use the mapping theorem fixing the form of the propagator in the infrared, 
e.g. in the Landau gauge, as 



DfM = Sab (fir 



p2 



n=0 



mi + ie 



- + 



1 



Ng 



(25) 



but this can be recomputed in any gauge by the classical equations with the mapping 
theorem. The next step is to use the approximation that holds in a strong coupling 
limit 



A; = aJ d 4 x'Df v {x - x')t{x') + O (-^j + 0(f ) 



(26) 



and we note that, in this approximation, the ghost field just decouples and becomes 
free and one finally has at the leading order 



Zo[j] = exp 



% - I d 4 x'd 4 x"j a ^(x')D^ u (x' - x")j bv {x") 



(27) 
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This functional describes free massive glueballs that are the proper states in the 
infrared limit. Yang-Mills theory is trivial in the limit of the coupling going to infinity 
and we expect the running coupling to go to zero lowering energies. Now, we can take 
a look at the ghost part of the action. We just note that, for this particular form of 
the propagator, inserting our approximation into the action produces an action for a 
free ghost field. Indeed, we will have 



c a d^c a + O 



Ng 



+ o(f) 



(28) 



So, a ghost propagator can be written down as 



p 2 + ie 



Ng 



(29) 



Our conclusion is that, in a strong coupling expansion 
decoupling solution. 



1/y/Ng, we get the so called 



4 QCD at infrared limit 



When use is made of the infrared fixed point result, QCD action can be written down 
quite easily. Indeed, we will have for the gluon field 



S gf = ^ / dVdV 



and for the quark fields 
Sg = J d 4 xq(x) 



r{x')Dt{x'-x")f\x") + 



1 



Ng 



+ o (f) 



(30) 



10 -m q - gY 



d A x'Df v (x - x')f\x' 



g 2 r 



M y / d 4 x'D;l(x - x') ^V)jYqV) + O + O (f) 



(31) 



q(x) 



We recognize here an explicit Yukawa interaction and a Nambu-Jona-Lasinio non- 
local term. Already at this stage we are able to recognize that NJL is the proper 
low-energy limit for QCD at zero temperature. Now we operate the Smilga's choice 
rfftlv = dabiVtw ~ VnVv/v 2 ) f° r the Landau gauge. We are left with the infrared limit 
QCD as one has using conservation of currents 



>9f 



X - J d 4 x'd 4 x" 



jZ(x')G(x'-x")r(x") + 



Ng 



+ o(f) 



(32) 
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and for the quark fields 

S q = J2 [ d A xq(x) i$-m q -g^Y I ' d 4 x'G(x - x')j^(x') (33) 
q J V J 

- sVy / d'x'G{x - x') E sV) yWM + (-^j + (f) 



q[x). 



This action can be the starting point for our analysis at finite temperature. But before 
doing this, we want to give explicitly the contributions from gluon resonances. In order 
to do this, we introduce the bosonic currents j°(x) = r}°^j{x) with the current j(x) 
that of the gluonic excitations after mapping. Using the relation rj^rj^ 0, = 3(N^ — 1) 
we get in the end 



S gS = \{Nl-l) j d'x'd 
and for the quark fields 

Sq = E j rf4 ^( X ) 



j{x')G{x' - x")j{x") + O 



Ng 



+ o (f) 



(34) 



a u 



A" 



i(p-m q - ffqfl 



<rx'G(x-x')j{x') 



(35) 



- sVy J d 4 x'G(x - x') J2 q'W)^q'(x') + O (-±=- 



q(x). 



Now, we recognize that the propagator is just a sum of Yukawa propagators weighted 
by exponential damping terms. So, we introduce the a field and truncate at the 
first excitation. This is a somewhat rough approximation but will be helpful in the 
following analysis. This means that we can write the bosonic currents contribution 
as coming from a boson field written down as a(x) = y^3(N^ — 1)/B J d 4 x'A(x — 
x')j{x'). So, the model we consider for our finite temperature analysis, directly derived 
from QCD, is 



d x 



1 

-777,5(7 



2„2 



(36) 



and for the quark fields 



S q = Y J d A xq(x) 



10 — rn q - 9\ 



(37) 



3Vy J d'x'G{x - *0 E^')y7^V) + O (-^L) + O (f) 



q{x) 
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Now, we are able to recover recover the non-local model of Weise et al. [12] directly 
from QCD (2Q(0) = G is the standard Nambu-Jona-Lasinio coupling) by setting 



G{p) 



Br, 



G 



n=0 



p 2 - (2n + l) 2 (ir/2K({)) 2 a + ie 2 



C(p) 



(3? 



with C(0) = 1 fixing in this way the value of G using the gluon propagator. 
We move to an Euclidean action and define the following fields 

<j) a (x) = (a(x),7r(x)) 

So, the bosonic action will be, after taking the expansion around the v.e.v. <ft a 

l^, 1 



S 



B 



d x 



-{d8af - -ml(5af 



+ S MF + S (2) + . 



1 



(39) 

= M), 

(40) 
(41) 
(42) 

m 2 + l/Gj (2n)^^ J p 2 + M 2 (p)' ^ 
The next-to-leading order term is given by (a correction to mass mo for the a field) 



gl M/ + M>(p)] + i 



being 

S M f/V 4 = -2NN f 
This holds together with the gap equations 

M(p) = m q + C(p)v 

and 



2 \G 



+ m ]v 



ANN 



f 



d*P M(p) 



S (2) = i 

2 



d A p 



[F + (p 2 )5a(p)5a(-p) + F_(p 2 )57r(p)57r(-j9 



(44) 



being 
^ ± (P 2 ) 



G 



4NN 



iS C( q )C( q+P ) ,«±tt TM y M i q + P) ^ (45) 



(2tt 



[q 2 + M 2 (q)][(q + p) 2 + M 2 (q+ p)] 



For the chiral condensate one has 



-ANN 



d A p 



M{p) 



p 2 + M 2 (p) p 2 + m 



(46) 



Till now there are two novelties really implied with respect to the work of Weise 
et al.[12\: The model is exactly obtained from QCD and the expression of the form 
factor C(p) is properly fixed through the exact gluon propagator at infrared. The next 



S 



step is to consider the case of finite temperature. This can be easily accomplished 
with the exchange 

being the sum over k that on Matsubara frequencies u k = 2kir/P for bosons and 
Uk = {2k + l)7r/(3 for fermions. So, we can write down the gap equations at finite 
temperature as 

M{cu k , p)=m q + C(cu k , p)v (48) 
ANN t o-X ^ f d3 P nr.. _n M(u k ,p) 



ml + l/G 



k=~oo 

while for the chiral condensate one has 



7g" i-J Wf c(ut,p) ^wrm^F) (49) 



m - -^iv- 1 I L.^L 3 - , m j. _, 1 • (50) 

? — „ 



(2VT) 3 L^+P 2 + M 2 ( Wfc ,p) W 2 +p2 + m 2 



Assuming the integral regularized by a cut-off A and noting that C(p) is practically 
1 in the low-energy range, we can prove the existence of a critical temperature where 
the chiral symmetry is restored. Setting v — into the gap equation we have to solve 

4NN f n , ^ f <$p s 1 

' ' V ' C 2 ( Wfc ,p)^— = 1. (51) 



' X 1 



mi + l/G" ^7 (2^)3" l — 
At small temperatures we are able to get the critical temperature 

7T 2 / „ 1 



^2 3 

T 2 « — 

c _2 



A2 - 7v7v7 l m « + c 



(52) 



This shows, starting directly from QCD, that a critical point does exist for this theory. 
We note that for Nf = 2 and T c = 170 MeV gives A = 769 MeV, perfectly consistent 
with NJL model. This expression is very similar to the one obtained in [22J. 

For aims of completeness, we give in figJH a comparison of our gluon propagator 
(the form factor) with that used in Weise et al. [12J and the exact one for the instanton 
liquid as given in [23] that inspired Weise et al. This shows that the approximation 
of a QCD vacuum as an instanton liquid is a very good one. 

I am much grateful to Marco Ruggieri without whose help I would not have 
obtained these results. 
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Figure 1: Comparison between the form factors with respect to the case of an instan- 
ton liquid. 
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